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Abstract: We provide an explicit construction for Gazeau-Klauder coherent states re- 
lated to non-Hermitian Hamiltonians with discrete bounded below and nondegenerate 
eigenspectrum. The underlying spacetime structure is taken to be of a noncommutative 
QJ , type with associated uncertainty relations implying minimal lengths. The uncertainty 

relations for the constructed states are shown to be saturated in a Hermitian as well 
as a non-Hermitian setting for a perturbed harmonic oscillator. The computed value 

T 1 

J> ' of the Mandel parameter dictates that the coherent wavepackets are assembled accord- 

ing to sub-Poissonian statistics. Fractional revival times, indicating the superposition of 
fNJ ■ classical-like sub-wave packets are clearly identified. 

cn 

i> 
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1. Introduction 

Noncommutative spacetime structures are suggested by gravitational stability [|J in almost 
all promising approaches to quantum gravity, such as string theory [||, [|, ||] or loop quantum 
gravity ||, ^| as well as black hole physics J?]]. Besides the numerous possible structures 
we focus here on a particularly interesting one giving rise to minimal measurable distance 
beyond which the entire concept of length becomes meaningless. Such type of cutoff in our 
possible knowledge of space results from generalized versions of Heisenberg's uncertainty 
relations, i.e. from modifying standard spacetime to certain noncommutative versions in a 
specific way. These type of spaces have attracted considerable attention in recent years in 
the mathematical and physical literature ||, |, |T|, [TT], HI H, HI H, HI, Hl> EDI H Hi- 



So far most of the attention has been paid to the algebras, but fairly little, in compari- 
son to standard systems, to the actual nature and properties of the states it acts on. Here we 
focus on the explicit construction of coherent states respecting these modified uncertainty 
relations. Generic expressions for coherent states were propose by Gazeau and Klauder 
(GZ) |^2|, By construction the states are expected to be stable when evolved in time, 
in the sense that they remain coherent during the evolution process. For a one dimensional 
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harmonic oscillator on a noncommutative space we present here a computation of these 
states in first order perturbation theory. When dealing with noncommutative spacetime 
structures of the aforementioned type one also encounters an addition complication due 
to the fact that the associated canonical variables are in general no longer Hermitian with 
regard to standard inner products [10|. Consequently Hamiltonians formulated in terms 



of these variables are also no longer Hermitian. We adopt here a recent approach to non- 
Hermitian systems [24, 25, 26, E?], 28] which renders them self-consistent and physically 



meaningful. 

A striking feature of the coherent states presented here is the well known fact that 
in a certain parameter regime the original wave packet can be fully reconstructed after a 
specific time and can be interpreted as a superposition of classical-like sub- wave packets. In 
the system considered here the existence of these structures is a signature of the spacetime 
deformation and disappears in a standard setting. This indicates the interesting possibility 
that the structures of noncommutative spaces could actually be probed experimentally. 

Many of the computations presented here have been attempted previously by Ghosh 
and Roy in [29], but almost all our results disagree with their findings and conclusions 



which we believe to be conceptually and compositionally incorrect as we will point in the 
course of our presentation. 

Our manuscript is organized as follows: In section 2 we assemble various generalities on 
GZ-coherent states and show how the construction needs to be altered for a non-Hermitian 
setting. We also set up our notation for a standard perturbative treatment. In section 3 we 
construct GZ-states and evaluates various expectation values for a Hermitian as well as non- 
Hermitian setting, which we use to test the Ehrenfest theorem and the equivalence principle. 
In section 4 we present our analysis for the revival times structure. Our conclusions are 
stated in section 5. 

2. Perturbative GZ-coherent states for non-Hermitian Hamiltonians 

We commence by establishing our notations by collecting some well known facts about 
GZ-coherent states and indicate the necessary modifications needed for a non-Hermitian 



setting. The GZ-coherent states [22, ^j] for a Hermitian Hamiltonian h with discrete 



bounded below and nondegenerate eigenspectrum are defined as a two parameter set 

|J ' 7 '^ = WU) £ Jn/2eXT *Z ijen) \<t> n ) , J€*Z,T€*. (2-1) 

\ I n=0 V Pn 

The states \<p n ) are the orthonormal eigenstates of h, that is h\<p n ) = hue n \(j) n ). The 
probability distribution and normalization constant are 

p n :=l[e k and AA 2 (J):=^— , (2.2) 

k=i k=o Pk 

respectively, where the latter results from the requirement (J, 7, (j) | J, 7, <f>) = 1. 

Here we will consider Hamiltonians H on noncommutative spaces which are non- 
Hermitian with regard to the standard inner product. The construction for the coherent 
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states is then easily adoptable when we assume that the Hamiltonian H is pseudo/quasi 
Hermitian, i.e. the non-Hermitian Hamiltonian H and the Hermitian Hamiltonian h are 
related by a similarity transformation h = rjHr]~ l , with r\^r\ being a positive definite oper- 
ator playing the role of the metric. The corresponding eigenstates |$) and \<f>) of H and h, 
respectively, are then simply related as |<5) = r?" 1 \(f>). As observables are expected to be 
Hermitian, we need to change the metric when computing expectation values for operators 



associated to the non-Hermitian system p4, 25, 26, E7L 28]. The same reasoning has to 



be adopted for the evaluation of expectation values with regard to the coherent states. 
Therefore the expectation value for a non-Hermitian operator O related to a Hermitian 
operator o by a similarity transformation o = rjOrj -1 is computed as 

(J, % ^\0\J^,^) v := (J, 7, *|t,V? I J, 7,$) = (J, 7, 0| o I J, 7, ^>. (2.3) 

Our notation is to be understood in the sense that in the state | J, 7, (f>) and | J, 7, $) we sum 
over the eigenstates of the Hermitian Hamiltonian h and non-Hermitian Hamiltonian H, 
respectively. These states are continuous in the two variables (J, 7), provide a resolution 
of the identity, are temporarily stable, in the sense that they remain coherent states under 
time evolution, and satisfy the action identity 

(J,7,$|fr|J,7,^ = (J,7,^|/i|J l7 ,^) = hu,J. (2.4) 

This identity ensures that (J, 7) are action angle variables |22|, p3| . 

The main purpose is here to consider a model on a noncommutative space with nontriv- 
ial commutation relations for their canonical variables giving rise to minimal uncertainties. 
It is then interesting to investigate how close the GK-states approach the minimum uncer- 
tainty product and eventually might even become squeezed states. Thus for a simultaneous 
measurement of two observables A and B in this system we need to evaluate the left and 
right hand side of the generalized version of Heisenberg's uncertainty relation 



AAAB > - 

~ 2 



{J, 1 ,<S>\[A,B]\J, 1 ,$) n . (2.5) 



The uncertainties are computed as AA = (J, 7, $| A 2 \ J, 7, $) — (J, 7, $| A \ J, 7, <I>)^ and 
analogously for AB. In order to test the quality of the coherent states, i.e. to see how 
closely they resemble classical mechanics, we may also test Ehrenfest's theorem for an 
operator A 

ihj t (J, 7 + tu,$\A\J,j + tu),<f>) v = (J,7 + tw,*| [A,H] \ J, 7 + toj, $) v . (2.6) 

We used in ( |2.6j ) the fact that the time evolution for the states | J, 7, <E>) is simply imple- 
mented as exp(—iHt/h) \ J, 7, <3?) = | J, 7 + toj, $) , see |2^, [23j . Specifying the operators A 
and B we will also test below the correspondence principle. 

Here we present a perturbative treatment of the above considerations around ho for a 
Hamiltonian decomposable as h = ho + hi, with /iq \ n) = e^ \n). According to standard 
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Rayleigh-Schrodinger perturbation theory the first order expansions of the eigenenergies 
and the eigenstates are 



e„ = 4 0) + (n| hi \n) + 0(r 2 ) and K) = | n ) + X/ 



hi \n) 
roj ToJ 



+ 0(r 2 ), (2.7) 



respectively. Wherever appropriate, we then simply use these expressions in (2.1) for our 
computations. 



3. GK-coherent states for the noncommutative harmonic oscillator 

We will now construct the GK-coherent states and various expectation values for the one 
dimensional harmonic oscillator 



10, 13| 



2m 2 



defined on the noncommutative space satisfying 

[X,P] =ih(l + fP 2 ) , X = (l + fp 2 )x, P = p. 



(3.1) 



(3.2) 



Here f := T/(mujh) has the dimension of an inverse squared momentum with r being 
dimensionless. We also provided in ([T^) a representation for the noncommutative variables 
in terms of the standard canonical variables x, p satisfying [x,p] = ih~. The ground state 
energy is conveniently shifted to allow for a factorization of the energy. The Hamiltonian 



in (|3.1| ) in terms of x, p differs from the one treated recently in [2S] as we take a different 



representation for X and P, which we believe to be incorrect in [29] even up to 0(r). 
The so-called Dyson map 77, whose adjoint action relates the non-Hermitian Hamiltonian 
in ( |3.1| ) to its isospectral Hermitian counterpart h, is easily found to be 77 = (1 + fp 2 )~ x l 2 . 
With the help of this expression we evaluate 

h = vHrj- 1 = ^ + —x 2 + ^r(j?x 2 + x 2 p 2 + 2xp 2 x) - hw ( - + -) + 0(t 2 ). (3.3) 
2m 2 4n \2 4/ 

Taking now ho to be the standard harmonic oscillator, the energy eigenvalues for H and h 
were computed to lowest order in perturbation theory [||, to 



E n = hwe Tl 



hum 



+ 0(t 2 



(3.4) 



According to ( [2.7|) we now calculate the first order expression for the wavefunctions to 

(3.5) 



71} 



— ^(n - 3)4 |n - 4) + —Jin + 1) 4 \n + 4) + C(r 2 ) 
lb lb 



where (x) n := T(x + n)/T{x) denotes the Pochhammer symbol. We stress that it is vital 
to include the second and third term in \<p n ) in order to achieve an accuracy of order r for 



expectation values. In |29], where a similar computation was attempted, these terms were 
incorrectly ignored. The expression for E n coincides with the one found in p9| for r — > 2A, 
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as this computation only involves (k\ h\ \n). Given e n as denned by by the relation (|3,4|), 
we compute the probability density and the expansions of its inverse 



-r n 



and 



1 

Pn 



1 



3 + n 



4(n- 1) 



TV. T + 0{T 



(3.6) 



We use the latter expression to evaluate the normalization constant in ( |2.2| 



N 2 {J) 



l- T j- j* +0(t") 



(3.7) 



We have now assembled all the necessary quantities to define the GZ-coherent states | J, 7, <p) 
in fl2,ip and are in the position to verify the validity of some of the crucial requirements 
on them, test their behaviour and compute expectation values. 

, ^?], 28 ] , in a non-Hermitian setting the observables are not 



As is well known [24, 25 



dictated by the Hamiltonian and it becomes a matter of choice to select them or equivalently 



the metric [24]. In fact, this is also true for a Hermitian system, where, however, the choice 
of the standard metric seems to be the most natural one. Here we are mainly interested in 
the Hamiltonian H of Q3.1| ) with X and P as observables, but it will also be instructive to 
consider first the Hermitian system described by h with x and p being the observables of 
choice. 



3.1 Obervables in the Hermitian system 



At first we consider the Hamiltonian h in (|3.3| ) as fundamental and treat the variables x 
and p as observables in that system. Expectation values are then most easily computed 
by taking the states |n) to be the normalized standard Fock space eigenstates of the har- 
monic oscillator with usual properties a) \n) = \/n + 1 \n + 1) and a \ n) = \Jn\n — 1). To 
first order in r, we then compute the expectation values of the creation and annihilation 
operators 



(J, 7, <P\a\J, 7, 0) = V^7e-*T [l-I(2 + J + 4i 7 (l + J))] + Jj 3 / 2 e 3 ^ + 0(r 2 ), (3. 
(J, 7, <P\ ^ I J, 7, 0) = yfj& [l-~(2 + J- 4*7(1 + J))] + \ J 3/ V 3 ^ + 0{t 2 ). (3. 



For the details of the computation we refer the reader to the appendix. In what follows 
we will often drop the explicit mentioning of the order in r, understanding that all our 
computations are carried out to first order. Using the fact that x = ^h/(2mw)(a) + a) 
and p = i\fhmuj/2(a^ — a), the expectation values 



(J,l,(f>\x\J,J, 
(J, 7, (j)\ p I J, 7, 



2Jh 



V2Jmujh 



cos 7 — T 

sin 7 + r 



cos 7 / sin 27 
7 sin 7 H h J sin 7 I 7 H — 



, (3.10) 



sin 7 / sin 27 
7 cos 7 h J cos 7 7 — 
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then follow trivially from (3.8) and ( |3.9[) . Expanding x 2 and p 2 in terms of a' and a, a 
similar, albeit more lengthy, computation yields 

h 



(J,7,0|x I J, 7, 



2mu 
hrnuo 



[l + 4 J cos 2 7 - r J (67 sin 2 7 + cos 2 7 + 2) (3.11) 

— r J 2 (47 sin 27 — cos 47 + 1)] , 

[l + 4 J sin 2 7 + t J (67 sin 2 7 + cos 2 7 - 2) (3.12) 
+r J 2 (47 sin 27 + cos 47 - 1)1 . 



These two expressions may be used to compute the expectation value of h, as defined in 
( |3,3j ), with regard to the GZ-coherent states. The remaining term in h only needs to be 
computed to zeroth order to achieve an overall accuracy of order r. We therefore calculate 

(J, 7, 0] p 2 x 2 + x 2 p 2 + 2xp 2 x I J, 7, 0) = h 2 (l + 4 J + 2 J 2 - 2 J 2 cos 4 7 ) + O(t). (3.13) 

Summing the contributions from fl3.10| ), (3.11) and ( p. 12 ), together with the required pre- 
factors to make up the Hamiltonian h, yields the action identity (2.4) as expected. We 
remark that this crucial identity was violated in [29]. 

Employing the above quantities we can also investigate how close the coherent states 
approach the minimum uncertainty product. Assembling the required expectation values 
we then obtain 

h 



Ax 2 = (J, 7, 0| x 2 I J, 7, 

Ap 2 = (J, 7, 0| p 2 I J, 7,, 
and therefore 



(J, 7, 0| x I J, 7, 
(J, 7, 0| p I J, 7,. 



2muj 
hmoj 



AxA P =^ + 0(t 2 ). 



[1 + r J(cos 2 7 - 27 sin2 7 )] ,(3.14) 
[1 - t J (cos 2 7 - 2 7 sin 2 7 )] , (3.15) 



(3.16) 



Thus the states \J, 7, 0) saturate the minimal uncertainty in a simultaneous measurement 
of x and p and therefore constitute squeezed states for all values of J and 7 up to first 
order in perturbation theory. 

Using (3.10) we also verify Ehrenfest's theorem ( |2.6| ) for the operators x 



(J, 7 + £w,0| [x,h] I J, 7 + tu,i 



ih — (J, 7 + tu, 0] x I J, 7 + tuj, 4 
at 

ih ituj o 2 

(J, 7 + iw, 0| — p H (px + x p + 2xpx) I J, 7 + to;, 

m 2 



(3.17) 



2Jw 



m 



sin 7 + r 



1 



- sin 7 + cos 7 I (J + 1)7 + - J sin 27 



and p 



d 



(J, 7 + tui, 0] p I J, 7 + tU, ( 



= ( J, 7 + iw, 0] [p, /i] I J, 7 + to; 
= ( J, 7 + tw, 0] — ihmu) 2 x — iruj(px 2 + x 2 p) \ J, 7 + iw, 0) 
= -iV2Jr^h 3/2 uj 3/2 ' 

For convenience we abbreviated here 7 := 7 + 



(3.18) 



cos 7 + — [(3 J + 2) cos 7 — 4( J + 1)7 sin 7 — 3 J cos 37] 
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3.2 Obervables in the non-Hermitian system 

As stated, the system we actually wish to investigate is described by the non-Hermitian 
Hamiltonian (|3.1| ) with a non-trivial commutation relation (|3.2| ) for its associated observ- 
ables X and P. In order to test the inequality (|2.5|) we need to compute 



(7,7,^1X17,7,$)^ 
(J, 7, $| X 2 | J, 7,3% 



2J/i r 



cos 7 + — sin7(Jsin27 — 27(1 + J)) 



(3.19) 



[l + 4J cos 2 7 + r[l + J(2 - 2 cos 2 7 - 67 sin 2 7 ) (3.20) 



2mw 

+ 2J sin 2 7 (sin 27-27)]] 



We note here that the actual computation has been carried out by translating first all 
quantities to a Hermitian setting and then following the same reasoning as in the previous 
subsection. Combining ( 3.19[) and ( |3.20| ) then yields 



AX 2 = (J,7,<J?|X 2 |J, 7 ,<D) -(J, 7 ,$|X|J,7,$) 



(3.21) 



h 



2mw 



[1 + r (1 + J(2 - 27 sin 2 7 - cos 27))] 



The computation for the expectation values of P is simpler, since the metric commutes 
with p, such that 



(J,7,$|P| J, 7, $>„ = (J, 7, 0b I J, 7, 0) , and (J, 7, $| P 2 | J, 7, *)„ = (J, 7, 01 P 2 1 4 7, 0) , 

(3.22) 

and therefore 

AP 2 = Ap 2 . (3.23) 



Expanding finally ( ggjj ) and flggg ), we obtain 



MAP = - r 
2 



l + ^(l + 4Jsin 2 7 )] = |(l + f(J, 7 ,cD|P 2 |J,7,cI>)). 



(3.24) 



This means that also in the non-Hermitian setting the minimal uncertainty product for the 
observables X and P, commuting as specified in fl3.2|) , is saturated. Thus to first order in 
perturbation theory also the GZ-coherent states |J, 7, <&) are squeezed states, remarkably 
this holds irrespective of the values for J and 7. 



Apparently this result was also obtained in 2£], but our disagreement with the results 
presented in there is at least fourfold. Firstly, the authors used the incorrect representation 
for the canonical variables X and P as mentioned earlier. Secondly the authors computed 
conceptually the wrong expectations values even when using their representation. Thirdly, 
the authors only take the first order in ( |3.5[) into account and therefore miss out various 
terms contributing to the first order calculation in r. Finally, even ignoring the previous 
three points and adopting all the wrong concepts used in ]29|| , we disagree on a purely 
computational level with many of the expressions presented in there. 
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Next we also verify Ehrenfest's theorem (^3) for the operators X 

ih ~k ^ 7 + $ ' x ' J ' 7 + tw ' = ^ 7 + $ ' ^ ' J ' 7 + tuJi *^ ' ^ 3 ' 25 ^ 



( J, 7 + tw, $| -(P + f P d ) I J, 7 + tu, *)„ 



2Jw 



sin 7 + r 



(J + 1)7 cos 7 + - sin 7(2 + J — 3Jcos 27) 



and the operator P 



ih— (J,7 + tu,$\P\J,'y + tu),$) v = (J,7 + tw,$| [P,F] |J,7 + iu;,$) 



-W2Jmh 3/2 uj 3/2 



cos 7 + — [(3 J + 2) cos 7 — 4( J + 1)7 sin 7 — 3 J cos 37] 



(3.26) 



Taking now for simplicity 7 = 0, differentiating fl3,25| ) once again and combining it with 
(3.26) we obtain the corresponding identity to Newton's equation of motion 



(J,tu>,$\ X \J,tu),$) 



-lu 2 { J, to;, $| X + — {3XP 2 + 3P 2 X + 2PXP) j J, tu, *)_ , (3.27) 



The relations fl3,25|) and fl3.26|) were not recovered in p^| , where the comparison between 
the left and right hand sides mismatched. Instead of ( p. 27 ) the authors proposed a "cor- 
respondence principle with twist". According to our argumentation this is incorrect and 
there is in fact no reason to assume the Newton's equation is simply the same as the one for 
the standard harmonic oscillator. The reason for the discrepancy are the aforementioned 
conceptual and computational mistakes in 2S]. 



4. Fractional revival structure 

Further insights into the comparison between the classical and quantum description can 
be obtained from the revival time for wave packets. For a general wave packet of the form 
ifi = c n<ft n sufficiently well localized, i.e. being governed by sub-Poissonian statistics, 
near a submode n = n with energy E n , it was argued in that beside the revival of the 
classical-like wave packet after the classical period T c \ = 2irh/ \E' n \ one may also encounter 
so-called fractional revivals. These partial revivals of the original wave packet occur at 
times p/q(T rev ), with coprime integers p, q and the revival time given by T rev = AttH/ \E'^\. 
Depending on the values of p and q one can interpret the emerging features as different 
types of superpositions of classical-like wave packets. 



For the case at hand we may follow [23| and expand the energy E n in the expression 



for the wave packets ( 2.1 



\J,u)t,4>) = } y Cn(J) exp(-itE n /h) \4> T: 



(4.1) 



n=0 
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with weig hting function c n (J) = J n / 2 /N{ J)y/p^, about fi := (n) = Jd\nN 2 {J) / dJ . To 
first order in perturbation theory we then easily compute 

(n) = J-T^J + ^pj+ 0(t 2 ), and <n 2 > = J + J 2 - r ( J + 3 J 2 + J 3 ) + 0(t 2 ), (4.2) 

2 \ -{ n ) 2 = J- t (J+J 2 )+0(t 2 ). (4.3) 



such that 



An = (n 

Consequently the Mandel parameter Q |31|] turns out to be negative 



«:-$-i~£ + °C)«>, 



(4.4) 



suggesting a sub-Poissonian statistics. This implies that we have a strong localization 
around n required for the possibility of the revival of the classical-like sub-wave packet. 
According to the above expressions we compute the classical period and the revival time 
to 



2vr r 



1 + 2J)tt, and T ri 



47T 



UJ U) LOT 

We now use these quantities to analyze the behaviour of the autocorrelation function 



(4.5) 



A(t) := | (J, 7,0 \J,l + tu, 



(J, 7, $ |J,7 + tw,$) 



(4.6) 



In order to find a set of meaningful values for our free parameters J, r and also to find 
an appropriate upper limit cutoff in the sum ( [4.1[) , let us first investigate the weighting 
function c n (J). 




5.0 n 



Figure 1: (a) Weighting function for r = 0.1 with (n) = 1.24,2.25,3.04 for J = 1.5,3,4.5, 
respectively and (b) r = 0.01 with (n) = 2.93, 5.76, 13.72 for J — 3,6, 15, respectively. 

For the chosen values we observer in figure [l] that the wave packets are well localized 
around n resulting from (|4.2|) , such that the prerequisite for the validity of the analysis in 



|3C] is given. Increasing the values of J for fixed r we observe negative values for \c(J)\ 



2 



for large values of n, which clearly indicates that our pertubative expressions are no longer 
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Figure 2: (a) Autocorrelation function as a function of time for J = 1.5, t = 0.1, U) = 0.5, h = 1, 
7 = 0, T c i = 10.05 and T ICV = 251.32; (b) Autocorrelation function as a function of time for J = 6, 
r = 0.01, ui = 0.5, h= 1, 7 = 0, T c] = 11.74 and T rov = 2513.27. 



valid in that regime. We also note that n ~ 50 will be a sufficiently good value to terminate 
the sum in the expression for the autocorrelation function (4.1) analyzed in figure 

In panel (a) of figure || we clearly observe local maxima at multiples of the classical 
period T c \. As explained in the first full reconstruction of the original wave packet is 
obtained at T rev /2 which is clearly visible in panel (a). The fractional revivals are better 
observed for smaller values of r as depicted in panel (b). In that scenario the classical 
periods are so small as compared to the revival time that they are no longer resolved. We 
clearly observe a number of fractional revivals, such as for instance T rev /4 corresponding 
to the superposition of two classical-like sub-wave packets and others as indicated in the 
figure. 

Notice that our expression and our analysis presented here differs once again from the 



one in [29], where for instance the mandatory revival at T rev /2 was not observed. 



5. Conclusions 



Our central results is the construction of explicit expressions for the GZ-coherent states 
for a non-Hermitian system on a noncommutative space leading to a generalized version of 
Heisenberg's uncertainty relation. We showed that these states are squeezed for all values 
of J and 7 as they saturate the minimal uncertainty. Crucially we established two of the 
nontrivial GZ-axioms are satisfied. First of all the states are shown to be temporarily 
stable, i.e. they remain coherent under time evolution, and second the states satisfy the 



action identity (2.4) allowing for a close relation to a classical description in terms of 
action angle variables. We also demonstrated that when using the appropriate metric 
Ehrenfest's theorem is satisfied for the observables X and P. The desired resemblance 
of the coherent states with a classical description was further underpinned by an analysis 
of the revival structure exhibiting the typical quasiclassical evolution of the original wave 
packet. It should be noted that in the considered case the wave packet revival time ( fO| ) 
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depends explicitly on the deformation parameter r, such that a possible measurement could 
distinguish between a noncommutative and a standard commutative space. For instance, in 



the order of femtoseconds half and quarter revivals have been observed experimentally [32] 
for molecular wave packets described by anharmonic oscillator potentials with eigenenergies 
similar to fl3.4j ). Our analysis holds to first order perturbation theory in r and of course it 
would be very interesting to extend this to higher order or eventually to the exact case. 

There are various other directions into which our analysis might be taken forward. For 
instance, different types of models might exhibit a varied behaviour. More challenging is a 
construction for such states in higher dimensions. A systematic comparison with different 
types of coherent states would be insightful, especially with rare constructions related to 
non-Hermitian Hamiltonians |3 
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A. Appendix 



We present here a sample computation in order to make our results transparent and re- 



producible. In addition, we highlight the differences compared to [29]. One of the simplest 
computation in this context is to evaluate the expectation value of the annihilation oper- 
ator (|3.8| ). The evaluation of expectation values for different types of operators is more 



involved, but goes along the same lines. Using the expression for (2A) we compute 



<J, 7 ^|a|J, 7 ^) = -L f; j(m+n)/2eX Pj!2j em ^(^|a|^). (A.l) 



n,m=0 



With the expansion of \(f) n ) to first order in r we obtain 



Vnd~m,n-i + ^ (V (n + l)iVn - y/ (n - 3) 4 Vn - 4) <5 mi „_ 5 (A. 2) 

+ (V ( n + ^hVn + 4 - \/{n- 3)4^) <W»+3, 
\/n8 m , n -l + -\/{n + l)(n + 2){n + 3)5 m , n+3 , 



such that 



(J, 7,0| a I J, 7,0) = ^2 



rt-l/2 ry(e„_i-e n ) 



J n+3 / 2 y/(n + l) 3 e i T( e "+ 3 - e ") 



n=0 



-+0(r 2 



(A.3) 

The last sum has been ignored in [29], but is an important contribution to order r. Using 
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e n -i 



■1 — tit and p n = p n -i&n the first sum in ( |A.3| ) is evaluated as 



—i~i jn— 1/2 g-iyriT 



-17 



n 



M 2 



E 



1 Pn-1 



l--(l + n + 4i 7 n) + 0(r 2 ),(A.4) 



-17 



00 

4/ 00 



E 

n.=0 



E 

t(1 + 4h)^ J 



Pn-l 



For the second sum in ( |A.3| ) we use e n+ 3 — e r 
such that it becomes 



E • 

n=0 



+ 0(r 2 ), 



1 - - (2 + J + 4*7(1 + J)) + 0(r 2 ). 



(A.5) 



(A.6) 



E 



jn+3/2 % /(n+l) 3e i7(3+3r(2+n)) 
4^ 2 PnV e «+3 e n+2en+l 



+ 0(r 2 



= 3 + 3r(2 + n) and p n+3 = p n e n+3 e n+2 e n+ i, 

-j3/2e^ ™ Jn r j3/2, 



4AA 2 ^o n 

n=0 "■ 



jn 



,3i7 



+ 0(r 2 ). 
(A.7) 

Collecting ( A.6 ) and ( A.7 ) we obtain ( |3.8[) . Similarly we compute the expectation values 
for x 2 , p 2 , x 2 p 2 etc by converting them first into expressions involving the a and and 
then using the arguments from above. For the noncommutative scenario we convert first 
to a setting involving Hermitian operators. For instance, we use 



(J, 7, *| X I J, 7, $>„ = (J, 7, <t>\ x + -(p 2 x + xp 2 ) I J, 7, 4>) + 0(t 2 
and compute the right hand side as explained above. 



(Ai 
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